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Abstract 

Two-dimensional random surfaces are studied numerically by the dynamical tri- 
angulation method. In order to generate various kinds of random surfaces, two 
higher derivative terms are added to the action. The phases of surfaces in the two- 
dimensional parameter space are classified into three states: flat, crumpled surface, 
and branched polymer. In addition, there exists a special point (pure gravity) corre- 
sponding to the universal fractal surface. A new probe to detect branched polymers 
is proposed, which makes use of the minbu(minimum neck baby universe) analy- 
sis. This method can clearly distinguish the branched polymer phase from another 
according to the sizes and arrangements of baby universes. The size distribution 
of baby universes changes drastically at the transition point between the branched 
polymer and other kind of surface. The phases of surfaces coupled with multi-Ising 
spins are studied in a similar manner. 
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1 Introduction 



Studies on two-dimensional surfaces are currently drawing much attention in various fields 
of science, such as physics, chemistry and biology. In physics, for example, the statisti- 
cal property of triangulated surfaces is under intensive investigations in the context of 
two-dimensional quantum gravity. For a surface without coupled matter (pure gravity), 
a Monte-Carlo simulation by a dynamical triangulation(DT) method shows the fractal 
property, as predicted by such analytical theories as the Liouville field theory fl],|],|3|], the 
matrix model 0] and string field theory ||. When matter, such as scalar fields, Ising 
spins or Potts models, is put on a surface it is expected that the surface becomes unsta- 
ble toward forming a one-dimensional structure, like a branched polymer ||[7|]. We are 
interested in generating and classifying various kinds of surfaces. It is easy to imagine 
that there exist at least three kinds of surfaces, ( flat, crumpled surface and branched 
polymer), and wonder whether any surface can be classified into one of these three types. 
We would also like to investigate the nature of the transition between them. Furthermore, 
if the fluctuation of a surface can be visualized, it will be very useful for understanding 
the properties of surfaces and mutual transitions. 

One of the important quantities which characterize a surface generated by the dynam- 
ical triangulation method is the discretized local scalar curvature; 



where qi (a coordination number) is the number of triangles sharing the vertex i, and 
a 2 is the area of an elementary triangle. As is well-known, the surface of pure gravity 
fluctuates strongly, containing many branches and spikes. A large positive curvature is 
localized at each tip of the spikes, and negative curvatures are condensed around the roots 
of the branches. Two additional terms are included in the action in order to control the 
local fluctuation of the surfaces. One is (3 J ^fgR 2 d 2 x\ the other is a J ^Jgg^ v d^Rd u Rd 2 x 
in continuous space, where a and f3 are coupling constants. The former affects on the 
local curvature of surfaces; i.e. for a large positive f3 it suppresses large fluctuations of 
R ||, and as (3 becomes negative the surfaces become crumpled or branched polymers, 
depending on the sign of a. 

Since the smallest neck size of a branch is limited by the lattice constant, it is significant 
to consider the minimum necks. Only for the minimum neck can we locate its position 
uniquely. By measuring the area distributions and connectivities of each baby universe 
separated by the minimum necks we can distinguish branched polymers from other kinds 
of surfaces. 

This paper is organized as follows: In the next section we briefly describe how the DT 
method is performed with two higher derivative terms in the action. In Section 3 we give 
a short review of the method of minimum neck baby universes (or 'minbu' M). In Section 
4 we propose three phases to classify random surfaces, and discuss the properties of each 
phase. In Section 5 the transition between a crumpled surface (or a fractal surface) and 
a branched polymer is discussed, together with numerical simulation of random surfaces 
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coupled with multi-Ising spins. In the last section we give a summary of our numerical 
results and a discussion. 



2 Dynamical Triangulation with Higher 
Derivative Terms 

A study of two-dimensional quantum gravity offers not only a simple model of Ein- 
stein gravity, but also a general framework for providing the universal property of two- 
dimensional surfaces. Here, we employ the dynamical triangulation method [lOj (DT), 
which is known to give the correct answer, expected from analytic theories. 

In DT, calculations of the partition function are performed by replacing the path 
integral over the metric to the sum over possible triangulations of the two-dimensional 
surfaces by equilateral triangles. From the triangulation condition and the topological- 
invariance the total number of i-simplices (iVj) is related as 

3N 2 = 2Ni, (2) 
N 2 -N 1 + N = X , (3) 

where x is the Euler number. 

In practice we proceed with DT for a surface with an arbitrary genus as follows: 

1. Initialization. 

Prepare equilateral 3-simplices(tetrahedrons), and put them together by gluing tri- 
angles face-to-face, forming a closed surface with a target genus. Then repeat the 
barycentric subdivision by choosing a triangle randomly until the area of the uni- 
verse is reaches to the target size. This manipulation is equivalent to the so-called 
(1,3) move. 

2. Thermalization of configurations. 

For changing the geometry so as to bring the surface into thermal equilibrium, pick 
up a pare of neighbouring triangles randomly and flip the link shared by them(that 
is the flip-flop or (2,2) move). This move conserves the area and topology of the 
surface, and is known to be ergodic, i.e. any two configurations of a fixed area and 
topology are connected by the finite sequence of the (2, 2) moves. The manifold 
conditions are always checked by requiring that any set of 2-simplexes having a 
0-simplex in common should constitute a combinatorial l-bah^S 11 ). We also restrict 
the geometry so as to not allow singular triangulations, such as self-energies and 
tadpoles, in the dual graph. 

For the case of pure gravity the action is independent of the geometry, and all graphs 
generated by the flip-flop moves are accepted with equal weight as those of a canonical 
ensemble. In this case no dimensional parameter enters, and the surface is expected to 
be fractal. 
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Here, we introduce two additional terms in the action in order to produce various types 
of surfaces. There is a correspondence between a continuous theory(characterized by the 
metric tensor g^) and a discretized theory(characterized by the coordination number qA 
through a relation 

a 2 N ° 

Vgd 2 x T X> ••• ■ (4) 

6 i=i 

Two new terms are expressed in the language of discretized theory by making use of 
eqs.(0) and (g) as 



and 



J ^g^d^Rd u Rd 2 x^ = -i= ^ J2 ( R i ~ Rj) 



v / 3\ 

= "4 (2Erf - E (6 - ft)(6 - ) , (6) 

va\ i a <i,j> I 

where (• • ■) means the ensemble average and indicates a nearest-neighbour pair of 
vertices. 

These two terms control the order of the vertices. The former favors a flat surface 
(i.e. qi = 6), and the latter creates a correlation of curvatures between neighbouring 
pairs of vertices producing such as crumpled surfaces or branched polymers for positive 
or negative parameter a, respectively. 

From eq.(|6|), since the / yfgg^d^RdvR term contains the R 2 term, we thus rearrange 
it into the J2<ij>(® ~ Qi){6 — Qj)/QiQj term in numerical simulations instead of the r.h.s. 
of eq. (]?]). Then, the Lattice action is 

q _n v (6 ~ q l ? v (6-gQ(6-%) 

<J — P Lattice <~ a Lattice ' V ' ) 

i 1i <i,j> QiQj 

where the script Lattice indicates a dimensionless coupling, 



a 



„ a Lattice: 
07T Z 



3a 2 J_ 

P . nKf^Lattice „ „ P Lattice) • 

47r 36 

Since these two higher derivative terms (eqs.(|5|) and ffl)) have dimensions of a~ 2 and 
a -4 , respectively, they are irrelevant in the continuum limit. The typical short- distance 
scale introduced by the R 2 term is ~ a//?, within which the surfaces are controlled to be 
smooth, and we need to set (3 > A in order to maintain an effect on surface significant 
over area A. 
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When both terms are included, it may be possible to construct an effective theory 
coupled with a scalar field by introducing an auxiliary field (x) and defining a theory for 
the case (3 > by the action 

/ d 2 i^{^g^d, X duX + ^X 2 ~ ixR), (9) 

and for (3 < 0, 

— ct 1 

d^Vai^g^d.xdvX + ^x 2 + xR)- (10) 

We then integrate out the auxiliary field formally, giving 

/ d^gR^nrR, (11) 

obtaining the effective action, 

S = J ^d 2 i{ag^d^Rd v R + I3R 2 ), for \a\ < (12) 

in a perturbative expansion using the parameters ^ . In order to guarantee the perturbative 
expansion of eq. ([□]), a and (3 must be restricted to \a\ <C \(3\ □. It is important to notice 
that the ratio of the couplings, -p-, can be fixed to be finite, while both parameters a 
and (3 are taken to have infinite simultaneity in order to keep the effects significant in 
the continuum limit. Since the signs of the kinetic terms of eqs.© and fllCf ) must be 
positive, it is required that a < 0. We again discuss the validity of the model of eq.(§) in 
subsection 4.1. 



3 Analysis of the baby universes 

In order to discuss the stability (or instability) of two-dimensional surfaces, we consider 
the so-called minimum- neck baby universes, as were mentioned in section 1. A minimum- 
neck baby universe is defined as a connected sub part of a surface whose area must be less 
than half of the total area, and the neck is constructed by three links which are closed, 
non-self intersecting. A mother universe is defined as a minimum-neck universe whose 
area is greater than half the total area. 

Suppose two kinds of surface configurations(in Fig. 1(a)): one is a non-branched sphere, 
whose volume is A; the other is branched sphere which has a baby universe (volume B). 
Once the latter configuration is favorable, all universes branch out boundlessly. Then, a 
bottle neck becomes a minimum-length neck, which leads to lattice-branched polymers. 
Instabilities make the surfaces branched, and simultaneously the mother universe disap- 
pears. Whether there exists a mother universe or not is one of the differences between 
branched polymers and other kinds of surfaces. 

*Due to the restriction of parameters(a, /?), our numerical results cannot include the so-called Feigin- 
Fuchs-type action, J d 2 ^^/g(ag^ v d^xd^X ~ ix^)- 
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It is easy to show that the relationship between of the string susceptibility to the 
branched polymer instability [jllj , for the partition function with the total area constrained 
to A given by 

Z(A) ~ k a A<-\ (13) 

For the case of a branched configuration (Fig. 1(a)), we can estimate a lower bound for 
this contribution, 



2A 



f 3 dBZ(B)Z(A - B)B(A - B) ~ k a A 2 ^ 3 . (14) 

"3 

Since Z[A] is the total number of distinct graphs with A triangles and S 2 topology, BZ[B] 
is the total number of distinct graphs with B triangles and D 2 topology for a boundary 
loop length of 3. 

The dominance of the branched polymer configuration corresponds to the case of 7 > 0. 
Thus, it is natural to recognize that a positive 7 corresponds to the case c > 1. So far, 
since there have been no exactly solvable models with c > 1, it is important to measure 
the actual minbu distributions for various kinds of DT surfaces in order to understand 
the dynamics of a 2d surface. 

In practice, it is easy to identify all of the minimum necks, and their connectivities on 
real dynamically triangulated graphs without self-energies and tadpoles in a dual graph. 
The procedures are as follows: 

(I) Pick up one link. There are then two vertices at both edges of this link; each 
vertex connects with many other links. Among them(links) a common vertex may 
be contained; unless these three links form a 2-simplex (i.e. a triangle) we call these 
closed links a minimum neck. A smaller part of the two universes separated by this 
minimum neck is a minbu. Repeat this process until all of the links are checked. 
This is applicable to the configurations of a nest of baby universes. 

(II) The connectivities of each baby universe are easily found. The key property is that 
a neck always separates universe into two parts for a spherical topology. First, 
pick up an arbitrary neck; then, two separated universes ({71, U2) can be identified. 
Second, choose a neck which belongs to universe {71 (or {72); then, Ul (or {72) can 
be decomposed into two sub parts: {711,{712(or {721, {722). Repeat until all of 
the necks have been checked. After this, it is straightforward to know the area 
distributions of all universes. 

The minbu analysis can also visualize a triangulated surface as a network of baby 
universes connected by minimum necks (Fig. 1(b)). In the next section we apply this 
analysis in order to classify the surfaces. 

A measurement of the minbu distributions provides the string susceptibility (7,,) with 



high accuracy [12]. The probability of finding a minbu with area B in a universe having 



a total area of A is given by 



3BZ[B](A-B)Z[A-B] 
Pa(B) = ^ (15) 
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(a) 




Figure 1: (a) Non branched surface with volume A and a branched surface which includes 
one baby universe with volume B. (b) A surface generated by the DT method visualized 
as a network of baby universes. One of the examples of the transition between fractal(or 
crumpled) surfaces to branched polymers is also visualized. 



In order to obtain the frequency of finding a minbu (Na(B)) with volume B for the 
case of pure gravity, we substitute the KPZ formula into eq.(^), giving 

^( fl )~^-'(ff"'{l-(f)p, (16) 

where j s = — ~. 

It is easy to extract the 7 S from numerical data when the KPZ formula is applicable. 
It should be remarked that in the case of c 7^ the finite size corrections for eq . (|16D 
become important. Even in the case of one Ising spin coupled to gravity (correspond to 
c = |), the finite size corrections are not negligible, and without a proper correction it 
is not possible to define the correct 7 a (= |) with the same accuracy as in the case of 
c = O0. 



4 Phase Diagram 

Let us classify typical states of surfaces generated by the two additional higher derivative 
terms in the action. We would like to propose that any surfaces are constituted by a 
branched polymer, a crumpled surface and a flat one. We can intuitively explain how these 
kinds of surfaces are created. The term J2i(Q~Qi) 2 with a positive coefficient favors $ = 6, 
i.e. a flat surface. On the other hands, the term X)t(f> — qi) 2 with a negative coefficient 
favors crumple surfaces rather than branched polymers. The term J2<ij> ^ 6 ~ g ^g 6 ~ 9 ^ with 
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Figure 2: Plot of Q 2 versus a and (3. 



a positive coefficient also makes surfaces crumple, because the opposite-sign curvatures of 
adjacent vertices become attractive. On the contrary the term J2<ij> ^~~^~ 2 ^ with a 
negative coefficient makes surface branched polymers, because same-sign curvatures are 
attractive. For the root (the loop of negative curvatures) and tip (a hemisphere of positive 
curvatures), those curvatures having the same sign need to condense to form a branched 
polymer. 

In order to classify the surface, we have observed Q2 and Q3 which are ensemble- 
averaged quantities (c/ eqs.(U) and (§)), defined on the lattice as follows: 

Q 2 = (J2^^)/N 2 , (17) 
T 9< 

Q 3 = (E (6 - qi)(6 - qj) )/N, (18) 

Figs.|2] and ^| clearly show that there are three plateaus, which indicate that the DT 
surfaces can be classified into three characteristic states (i.e. crumpled, flat and branched 
polymer) depending on the two parameters which we have employed. 

We have also measured the size-distribution and arrangements of all minimum-neck 
baby universes on these three plateaus. The features of these three surfaces are itemized 
as follows: 
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Figure 3: Plot of Q 3 versus a and (3. 

1. Crumpled Surface (a > 0, (3 < a) : 

This surface looks jagged. The value of Q2 becomes much larger than in the fiat 
case, and that of Q3 becomes negative. A mother universe exists whose volume 
reaches the order of 60% or 70% of the total area, although the growth of large 
branches is highly suppressed(see Fig.|5|). All baby universes have a small volume 
(~ 0(1)). The number of minimum necks is about 8 ~ 9% of the total number of 
triangles. 

2. Flat Surface {(3 > \a\) : 

Almost all fluctuations are suppressed (i.e. Q2 ~ Q3 ~ 0), and the fractal dimension 
becomes about 2.0. 

3. Branched Polymer (a < 0,(3 < \a\) : 

One of the typical branched polymers is illustrated by the r.h.s. of Fig.|l|(6). The 
mother universe completely disappears. All baby universes are very small (~ 0(1)). 
The number of minimum necks is about 30 ~ 40% of the total number of triangles. 

At the N2 —>■ 00 limit of the branched polymer its surface can no longer be regarded as 
being a smooth continuum medium. In order to check this, we measured the resistivities 
(correspond to the complex structures defined on DT surfaces) of branched polymers 
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Figure 4: Normalized histogram of the resistivities @ for the case of branched-like sur- 
faces: a = —10,(3 = 0. 

based on ref. 0. If there exists a smooth continuum limit of this surface, the peaks of 
resistivity become narrower as the size increases. Fig.[| is a plot of a normalized histogram 
of the resistivities for branched polymers. As has been expected above, Fig.f| shows that 
there is no sharpening of the peaks as the size increases. That indicates that a branched 
polymer has no continuum limit. 

Furthermore we point out the characteristic property based on an intrinsic geometrical 
point of view of these three states, that is, the distributions of the number of boundaries 
(branches) versus the geodesic distances (D). A definition of the boundaries (branches) 
is given in the next subsection (4.2). Fig.^] supports the intuitive pictures that we men- 
tioned above. Especially, the dashed-dotted line (correspond to the crumpled space-time) 
indicates that the large branch structures are highly suppressed. This tendency coincides 
with that of of three-dimensional simplicial gravity [jT4f| . 



4.1 Flat states (i? 2 -gravity) 

Here, we consider in detail the almost flat states of DT surfaces. At first, we investigate 
the minbu distributions for a large positive (3 (a = 0) ]nj. In ref. || the partition function 
and the string susceptibility have been estimated to be the following asymptotic form for 
4 — > and A -> oo: 

Z KN ~ A*"- 3 exp(-64vr 2 (l - hf ■ ^)e^ A , (19) 
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Figure 5: Distributions of the number of boundaries (branches) versus the geodesic dis- 
tance for various states with double-log scales. The solid line represents branched poly- 
mer states, the dotted line represents the case of pure gravity (fractal states) and the 
dashed-dotted line represents crumpled states. 



where 

7™(^2 + l^ X + ^f. (20) 

Here, x is the Euler number of this surface, and d is the central charge of the matter 
fields. For the case of d = and x = 2(sphere : S 2 ), becomes zero. In Fig.^ we 
plot and compare the minbu distributions (Na{B)) versus B(l - §) with the double-log 
scale for ^Lattice = 100 and 200 with a total number of triangles 5000. Since the R 2 
term affects only on the short- distance property of surfaces, it is expected that eq.(|2HD 



is adaptable for a relatively small minbu region. Figjg shows a good agreement with 
numerical data(open circle and square), and an asymptotic formula (dashed-line) which 
is obtained by substituting eq.([T9|) into eq.(|15"D within the small B regions I. We note that 
the data points approach the line predicted for the pure-gravity case (7 S = — ~) in large 
minbu regions. This indicates that the R 2 term is irrelevant with respect to long-range 
structures. 

When we consider a part of a flat parameter region, < a Lattice -C ^Lattice, the model 



' Since this asymptotic form is only assumed to be valid for B ^> 1, we should consider the simplest 
type of correction for B i.e. £? 7 ~ 2 — * £? 7 ~ 2 (1 + % + ■■■) [ p^[ . Two dashed-lines in Fig.^| are as follows: 



B d 

1)) = o - zaiitf(i - 

where Co and C\ are fit parameters 



\n(N A (B)) =C Q ~2 \n{B(l - -)} + — ^ , (21) 
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Figure 6: Minbu area distributions (Na{B)) for R 2 — gravity with topology S 2 with the 
double-log scales. The total number of triangles is 5,000. B represents the area of a 
baby universe. The solid line represents the minbu area distributions for the pure-gravity 
case (eq. (|i~6l) ). The two dashed lines show the fit curve predicted by R 2 — gravity |§ with 
pLattice = 100 and 200. 



of eq.(|9[) corresponds to our numerical simulations. The string susceptibility (jg ) of the 
model @ was also obtained by ref. || in the asymptotic form for ^ — > and A — > oo, 



KN = 2 _ 24 + , J_ riM±v\ X_ d_ £ /24 + , A 

,s I 12 32vr a V r] la. 87ra_r? a y V P 



where rj = ■ ^ — d and a_ = ^(24 + r/ — Jr](24 + 77)). It is remarkable that the 

of eq.fl22D does not depend upon a, and is the same as 7 S of eq. (|20|) when d — 0. The 

plateau(corresponds to [ar| < /3) observed in Figs.|2] and ^| indicates the a-independence of 



As reported in ref. [|T^, T7\, the transition from the fractal configuration to the flat 
configuration was a cross-over. It can be realized by varying f3 with a = in our case. In 
Fig.|7| we plot Log(Q2) versus Log(/3 Lattice), comparing the semiclassical calculation (see 
appendix A) and numerical results with the total number of triangles (20,000). A semi- 
classical analysis shows the excellent agreement with the numerical data at /3 ~ and 
large j3. 



4.2 Fractal structures at a ~ (3 ~ 

The most significant transition point is at (a ~ (3 ~ 0). Here, no scale parameters exist, 
and the surfaces are expected to become a fractal. We now define the intrinsic geometry 
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L %o( lattice) 

Figure 7: Plot of Log w (Q2) versus Log w (j3 Lattice), comparing the semi-classical calcula- 
tion and numerical results with a total number of triangles of 20, 000. The filled circles 
represent the numerical data with a Lattice — and ^Lattice — 0.01 ~ 300; the thin line is 
a guide. The thick line represents a semi-classical analysis. (See Appendix A). 



using the concept of a geodesic distance between two triangles on a DT surface. Suppose 
that a disk which is covered within D steps from some reference triangle. Because of the 
branching of the surface, a disc is not always simply-connected, and there usually appear 
many boundaries which consist of S 1 loops in this disc. In order to show the fractal 
properties of the DT surface, the loop-length distribution function (p(L,D)) is measured 
by counting the number of boundary loops with the length L which make boundaries of 
the area covered within D steps. In ref. @] it has been predicted that the loop-length 
distribution is a function of the scaling variable, x = l/d 2 , in the continuum limit for 
pure gravity, where I is the loop length and d is a distance defined on the continuum 
surfaces. For simplicity we do not distinguish continuous variables (I, d) from lattice 
variables (L, D). In Fig.[8]the results of our simulation for various distances are compared 
to the theory for a surface with the size of 20, 000 triangles. The distributions with 
different distances show an excellent agreement with each other and the shapes of the 
numerical data and the theoretical curve are quite similar. At the same time, artifacts 
of a lattice at small x and the finite-size effects at large x are observed. We have also 



investigated for several sizes of surfaces, i.e. 100,000 and 400,000 triangles [16], and the 



same quantities have been given(but low statistics). These excellent agreements with the 
numerical results and the analytical approach make it clear that the DT surface becomes 
fractal in the sense that sections of the surface at different distances from a given point look 
exactly the same after a proper rescaling of the loop lengths. The loop-length distribution 
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Figure 8: Loop-length distributions with the double-log scales. The total number of tri- 
angles is 20,000, and is averaged over 1000 configurations, x = L/D 2 is a scaling variable, 
where D represents the geodesic distance (it was measured at steps D = 15, 20, 25, 30) 
and L represents each loop length at step D. The small circles, quadrangles and triangles 
indicate the results of numerical simulations; the solid line indicates the theoretical curve 
predicted by string field theory ||. 



function comprises two types of loops corresponding to the exponent: baby loops and a 
mother loop. The baby loops dominate the loop length, and give the fractal dimension 
4; however, their averaged length is on the order of the cut-off length, and they are non- 
universal, and it cannot be considered as the proper fractal dimension of this surface. 
There always exists one mother loop with length proportional to D 2 , which has a fractal 



dimension of 3. In our earlier simulations fLofl , we found that the naive fractal dimension 
reached about 3.5 with 400, 000 triangles. 



5 Transition to a Branched Polymer 

The a term (a > 0) favors the attraction of vertices with the opposite-sign curvature. If 
this coefficient still becomes larger, this effect exceeds the R 2 term effect, and the surface 
crumples. On the other hand, a negative coefficient (a < 0) makes a transition to the 
branched polymer. In Fig.|l|(6) we illustrate one such kind of transition. 

The important point to note is the volume distributions of a mother universe and a 
baby one. For fractal surfaces(corresponding to the pure-gravity case), a mother universe 
has a large volume i.e. ~ half of the total volume( 0(N 2 )), and the sizes of the other baby 
universes are very small (0(1)). On the other hand, for branched polymer configurations 
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-60.0 -40.0 -20.0 0.0 20.0 40.0 60.0 
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Figure 9: RV 2 plotted versus OLLattice- The total number of triangles is 2,000, and (3 = 0. 
The statistical errors were estimated to be about 10%. 

all of the baby universes have small sizes and no mother universe exists. We therefore 
propose a new order parameter (RVi^) for the transition to the branched polymer, 

x-,T t V max ^second /<->o\ 

RVi = - — or RV 2 = — , (23) 

''total ' max 

where V max represents the maximum area among all of the universes in this surface, and 
y second represents a secondary size universe. Indeed, we examined that RV\^ works well 
as an indicator of such kinds of transitions in numerical simulations. Fig.|| shows the RV 2 
versus otLattice with the range: [—50, 50]. A clear transition can be seen from the branched 
polymer to crumpled or fractal surfaces. 



In FigsJlQ and 11 we also show RV\, as well as the total number of minimum necks 



versus a Lattice, respectively. These data (Fig.p],pT0| and [TTJ) suggest that the DT surfaces 



become a branched polymer for a < —10. As has been mentioned, our numerical systems 
correspond to the effective theory coupled with a scalar field; also a must be negative (i.e. 
ot Lattice must be positive). Then, the branched polymer stands for the appearance of some 
instabilities of the system. In our results if a Lattice < means the negative coefficient of 
the kinetic term (^|), and the stability of the surfaces certainly disappear. 

5.1 Matter (c > 1) coupled with gravity 

There are many theoretical indications that the surfaces degenerate into branched-polymer 
configurations for c > 1. We investigated the case of multi-Ising spins coupled to gravity 
using the minbu analysis. Each Ising spin lives on a triangles, and is thermalized by using 



cluster updates |L9 
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Figure 10: RV\ plotted versus oiLauice- The other conditions are the same as in Fig.9. 
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Figure 11: Total number of minimum necks plotted versus a Lattice- The numerical data 
are indicated by filled circles. The other conditions are the same as in Fig.9. 
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Figure 12: i?V 2 plotted versus the coupling constant of the Ising spins (Rising)- The 
open diamond indicates one Ising spin(c = 0.5), the close diamond indicates two Ising 
spins(c = 1), the open triangle indicates 4-Ising spins(c = 2), the close triangle indicates 
6-Ising spins(c = 3), the open circle indicates 8-Ising spins(c = 4) and the close circle 
indicates 12-Ising spins(c = 6). The total number of triangles is 2,000, and averaged by 
500 ensembles for each data. 
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Figure 13: RVi (V tota i = 2000) plotted versus f3 Ising ] the other conditions are the same as 
in Fig.12. 
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Figure 14: Total number of minimum necks plotted versus flising, the other conditions are 
the same as in Fig. 12. 



It is rather difficult to obtain clear signals for such transitions when n-Ising spins are 
coupled to gravity. In Figs.[l2|, [13] and [L4] we plot -RV^i and the averaged total number 
of minimum necks of the surface as a function of the coupling strength of the Ising spins, 
respectively. All of these observables slowly change with the central charge; however, it 
is certain that the effect of matter makes the surfaces branched for a sufficiently large c. 

It seems natural that the Ising spins living on a small-size baby universes, or branched 
polymers, suffer considerably from finite size effects. It is thus expected that such spins 
decouple with the background geometries. This is the main reason that the transition is 
not very clear. 



6 Summary and Discussion 

We have classified the two-dimensional random surfaces generated by dynamical trian- 
gulation into three states (crumpled surface, flat surface and branched polymer) varying 
the strength of two interaction terms. The three states are characterized by local observ- 
ables, such as average square curvatures (Q 2 ) and the average square-curvature derivative 
(Q3). The R 2 — gravity in the strong coupling limit 4 — > 00, gives the flat-surface, which 
was established numerically. We could also explain the cross-over phenomenon between 
the fractal and flat states within the perturbative colculations. Furthermore, we exam- 
ined why the surfaces become crumple rather than a branched polymer in the region of 
ot Lattice > 0, P Lattice < & Lattice- On the other hand we observed that the surface becomes 
branched polymers with a Lattice < 0, /5 'Lattice < |ct ■ Lattice], and that the measurements of 
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its resistivities clarify that this surface does not have a smooth continuum limit. The 
minbu analysis showed the nested structure of random surfaces directly; it also gave one 
important clue for distinguishing branched polymers from other kinds of surfaces. As 
to whether a mother universe exists or not, we quantitatively propose an observable, 
RVi = Yr 1 ^ or RV2 = v ^ cond . In fact, the RV\ 2 works well in our simulations. The good 

Vtotal Vmax ' 

fractal structure which is known in the case of pure gravity (a ~ (3 ~ 0) disappears in 
each of three states; on crumpled surfaces large branches can not grow enough, and on 
branched polymers a mother universe disappears, while branching into many small baby 
universes, which means a disappearance of the mother loop in the distribution of the loop 
length. We also applied the same analysis to the case of multi-Ising spins coupled to 
gravity corresponding to c = 1,2,4,6,8, 12. The RVx^ smoothly varies, and the number 
of minimum necks slowly grows with the central charges. 

The kinetic terms of the auxiliary field of eqs.(^) and (|T0D have negative signs for 
(^Lattice < within which we observed the branched polymer configurations in our numer- 
ical simulation. Since we do not know a connection between the effective theory coupled 
to a scalar field and the numerical simulation clearly, the central charge (c) which cor- 
responds to the DT surfaces cannot be estimated precisely. However, we can be fairly 
certain that these branched polymers are related to the c > 1 instabilities of the sur- 
faces. We therefore need to inquire further into the intimate correspondences between the 
numerical simulations and continuum theory. 
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Appendix 



A Semi-classical calculation of i? 2 -gravity 

In this appendix we consider a semiclassical treatment based on a ref. |2(| of i? 2 -gravity in 
order to account for the cross-over transition between a flat surface and a fractal surface. 
The partition function for surfaces with a fixed area of A is defined by 

z w = IiMk e ~ SoU ~ SJXI <J dx ^- A )' (A1) 

where Sa[g] = J d 2 x^fg{\ + PR 2 ), and 5 m [X] is the matter action given by S m [X] = 
f d 2 x y fgg» v d lx X a d u X a : a = 1 ~ d. 
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Taking the conformal gauge = g^e^, and estimating the FP-determinant and the 
integral of matter fields, we have 

Z (A) = j ^ -S G[ ge*]-^f d > xV g(- < t>A 1< t> + 2R<p) § ( f faTfr* _ A ? (A3) 

J Vol. (C. G.J \J J 

where A± is the Laplacian of S 2 with radius 1. Then, integrating the constant mode of 
(j), which is only used to estimate the 5-function, as O — hi A — lnG[0], where G[4>] = 
J d 2 x^e^, we obtain 

Z{A,(3) = A d -^e~ XA J V § fe- § M, (A4) 
where the action for oscillating parts, = — (f> , is now 

S[$] = /rf 2 ^(^(2-A 1 0)V^-^Mi0-^lnG[0]). (A5) 

Here we have written 7 = Now, we make a weak field expansion, and take only 

the quadratic terms with respect to the field, 0. The functional integral is carried out by 
expanding in terms of spherical-harmonic functions. Finally, we obtain the expression 



2 L 

„, . , rf-20 % a 64^/3 -, r , s 21+1 , , 

Z(A,/3) = A— e- M — rm)—, (A6) 

l>2 

where 

., = (^ (i(i+1) _ 2r+ ^w±i)^). 

The zero-frequency modes with I = 1 correspond to the conformal Killing modes, and 
have been omitted in the functional integral. 
For a large (3, the determinant part is 

Ef> 2 (2Z + l)ln^~ ( 2 L(L + 2)-| + 34)lnL 2 -2L 2 

-2(L(L + 2) -3)lnA + (L(L + 2) -3)ln(167r 2 /3). 1 ' ' 

The ultra-violet cut-off of an angular momentum (L) will be fixed by L(L + 2) = AA 2 . 
We obtain 

L /in A \ 

52(21 + 1) lnw, ~ k + ki(/3)^ + — + — - In A - 3 ln(167r 2 /3), (A8) 
i>2 V 3 87r P7/ 

where k and are divergent when the cut-off (A) goes to infinity. 
For small f3, we have 



V(2/ + l)ln^ ~ K > + K ' 1 (p)A + l\nA + 3\n(-?-) . (A9) 
i>2 3 V2tt/ 



19 



In order to compare the semi-classical calculation with the numerical results, we should 
make a correspondence of the cut-off (A) with N2 as L(L + 2) = AA 2 = kN 2 , where k is a 
constant which should be determined based on the numerical data. Using the asymptotic 
formulas eqs.(A8) and (A9), we obtain (/ cPxR 2 } = — 81 af, Z for large ^Lattice, 



a 2 



yv5flw«)=«^ ar- + s^ +! ^ i ■ (A10) 

\ J / {^PLattice OiV 2 Phatticel J 

and for small ^Lattice 

a 2 (/ y/g&fx) = 12nN 2l K + i-j . (All) 

We are able to determine the value of k from a data point, for example (Q 2, P Lattice) = 
(0.0449 ±0.0005, 200). From eqs.(|),(|),(|17D and (A10) we obtain 

K — 0.499 i 0.005, (A12) 



which is consistent with ref. 16 



From eq.(A6) we can calculate (/ cPx^/g R 2 ) without using the asymptotic formulas 
given above. We thus have 

(J dQR 2 ) 4 1 ^ 2Z + 1 ^(Z(Z + 1) - 2) 

48n 2 NA 2 3N 2 N ~^ 2 (/(/ + 1) _ 2) + M ' 1 ' ' 

which can well reproduce the data (Fig.[?p in both ranges of small and large (3, when we 
choose the parameter k to be 0.5. 

The values of ^string given through semi-classical calculation, 

f for small /3 

lstrm9 = \ d -^+ { Ww ^ large/?, (A13) 

is exact in a large (-d) limit. If one wants to obtain correct values for all d, a full quantum 
treatment is necessary. 
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